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Abstract 

We show that one-dimensional quantum systems with gapless degrees of 
freedom and open boundary conditions form a new universality class of quan- 
tum critical behavior, which we propose to call "bounded Luttinger liquids" . 
They share the following properties with ordinary (periodic) Luttinger liq- 
uids: absence of fermionic quasi-particle excitations, charge-spin separation, 
anomalous power-law correlations with exponents whose scaling relations are 
parametrized by a single coupling constant per degree of freedom, K v . The 
values of K v are independent of boundary conditions, but the representation 
of the critical exponents in terms of these K v depends on boundary con- 
ditions. We illustrate these scaling relations by exploring general rules for 
boundary critical exponents derived earlier using the Bethe Ansatz solution 
of the ID Hubbard model together with boundary conformal field theory, and 
the theory of Luttinger liquids in finite-size systems. We apply this theory 
to the photoemission properties of the organic conductors (TMTSF)2X , and 
discuss to what extent the assumption of finite strands with open boundaries 
at the sample surface can reconcile the experimental results with independent 
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information on the Luttinger liquid state in these materials. 
PACS numbers: 05.70.Jk, 68.35.Rh, 71.10.Pm, 79.60.-i 
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I. MOTIVATION 



Phase transitions take place in a different way on surfaces and in the bulk of a sample^. 
Order parameters, critical temperatures, critical exponents, and their scaling relations may 
be different, i.e. entirely new universality classes may be realized. Moreover, one may observe 
new phenomena due to anisotropy and the breaking of translational invariance caused by 
the boundary, like oscillations in correlation functions which, in the bulk, are monotonous, 
coordinate dependences, and in particular Friedel oscillations, in local quantities, etc. 

Here, we discuss such boundary critical phenomena in one-dimensional (ID) strongly 
correlated electron systems which possess, in the bulk, a quantum critical point at zero 
temperature. Conformal invariancel!, a consequence of the combined Lorentz and scale 
invariance at a critical point, allows an exact determination of all critical exponents in 
1+1D. The conformal field theory is parametrized by a unique constant - the conformal 
anomaly or the central charge c of the corresponding Virasoro algebra!. Strongly correlated 
electron systems of the kind we are interested in here, have c = 1, and their critical exponents 
continuously depend on the coupling of the fields. They can be calculated exactly from the 
low-energy excitations of the underlying HamiltoniansBU. The excitation spectra, in turn, 
can be obtained by a variety of methods, notably Bethe Ansatz for integrable models, and 
exact numerical diagonalization on small lattices quite generally. 

Conformally invariant, ID strongly correlated electron systems are metals, and are also 
described as Luttinger liquidsHi. In that perspective, the main focus is on their non-Fermi 
liquid properties which are embodied in their critical exponents (anomalous dimensions) and 
different velocities for excitations in different conformally invariant sectors (for two sectors 
only, charge and spin, this leads to charge-spin separation). In Luttinger liquid theory, 
the scaling relations between the different critical exponents are parametrized by a single 
renormalized coupling constant per degree of freedom K v (here v — p, charge, and v = a, 
spin), playing the role of the Landau parameters familiar from Fermi liquid theory. The 
"Luttinger liquid universality class" then corresponds to unique dependences of all critical 
exponents on the K v . The connection between both methods is well established in bulk 
systemsili'i. 

Recently, both approaches have been extended to boundary critical phenomena in ID 
fermion systems with open boundary conditions. Boundary conformal field theory! has 
been used to derive the boundary critical exponents of integrable ID electron systems in 
terms of the dressed charge matrices of these mo (Following the early work of 

Gaudinlll, Bethe Ansatz solutions have been produced, and integrability proven, for certain 
ID electron models with boundaries^iB.) Also Luttinger liquid theory has been formulated 
for systems with open boundaries@~0. While Wang et a/.0 gave the rules for evaluating 
critical exponents, explicit expressions for the exponents of specific correlation functions, 
relevant values for specific models such as the ID Hubbard model, and their relation to 
those derived within Luttinger liquid theory, are still missing. An exception are Friedel 
oscillations in the ID Hubbard model which have been studied by Bediirftig et a/.Jll. Here, 
and for spinless fermionsti!, accurate DMRG calculations are in impressive agreement with 
predictions from conformal field theory0. Related problems have also been studied in the 
context of quantum spin chains^. 

Here, we provide explicit expressions for the critical exponents of a variety of correlation 
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functions, and explicitly connect the predictions of conformal field theory in terms of the 
dressed charge matrices of integrable models to standard Luttinger liquid notation. We 
discuss the consequences of such mappings for nonintegrable systems and for experiments. 
We discuss the consequences of broken translational invariance in open systems. To be 
specific, we discuss the ID Hubbard model throughout this paper, except when stated 
otherwise. The results can be carried over to other integrable models without difficulty 

Quite generally, we are interested in how quantum critical ID strongly correlated elec- 
tron systems fit into the general framework of surface critical phenomena, and of different 
universality classes in the bulk and at surfaces. We find that the boundary critical exponents 
of gapless (conformally invariant) ID electron systems define a new universality class which 
we propose to call bounded Luttinger liquids, and which affect almost all physical quantities. 

We are motivated by several recent developments, both theoretical and experimental. 
Recent theoretical work on boundary effects in photoemission emphasizes the need for exact 
information on the spectral function of ID Hubbard models with open boundaries. Ap- 
proximate results would indicate an enhancement of spectral weight in a large energy range 
around the chemical potential^. On the other hand, unexplained photoemission experiments 
on the quasi- ID Bechgaard salt^HJH have been tentatively associated with the possible influ- 
ence of boundary effectsBil. Finally, boundary critical exponents are relevant in mesoscopic 
quantum wires and quantum Hall edge states and carbon nanotubes&llzl, and the Luttinger 
liquid language has been preferred in all these articles. 

As we shall see, integrable models can be used to illustrate relations between bulk and 
boundary critical exponents which are valid beyond the realm of integrability, in non- 
integrable Hamiltonians and even in experimental Luttinger liquids, provided they exist. 
Remarkably, they allow to propose boundary effects as a possible resolution of the qualita- 
tive discrepancy of the photoemission spectra of the ID Bechgaard salts, and the information 
provided by almost all other experiments on these materials. 

The reminder of this paper is organized as follows. In Section [TT] we discuss the deter- 
mination of critical exponents of ID correlated electron systems with periodic and open 
boundary conditions. We show which quantities determining the critical exponents are in- 
dependent of boundary conditions, and which do depend on them. We do this for the ID 
Hubbard model where we discuss the results from the Bethe Ansatz solution and use bound- 
ary conformal field theory, and for Luttinger liquids. For Luttinger liquids, we rewrite all 
critical exponents in a form similar to conformal field theory, i.e. in terms of the coupling 
constants K u , and the quantum numbers which an operator inserts into a system. The rules 
derived here are used in Section [Hj to give explicit expressions for a variety of correlation 
functions of bounded Luttinger liquids, and to illustrate the new scaling relations found here. 
Section [TV] discusses the application of boundary critical phenomena to problems posed by 
photoemission experiments on a class on ID organic conductors, and Section |V| summarizes 
the open questions remaining. 



II. MODELS 

While our results are valid generally for interacting ID electron systems, we shall concen- 
trate our discussion on the ID Hubbard model as a prototypical example of an integrable (or 
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Bethe ansatz solvable, used synonymously) system, and on the ID Luttinger model. Here, 
we briefly describe how their critical properties are calculated, and how the solutions with 
periodic and open boundary conditions are related. 



A. The ID Hubbard model 



The ID Hubbard model is defined by the Hamiltonian 

N-l N N r N 



ff = -E E (clc i+la + H.c.j + uJ2 n n n n -/iEE n » - o 5-^ T ~ Uil ^ ^ 2 -^ 

i=l cr=± i=l i=\ a=± j=l 

where C ia (Cj a ) is the electron annihilation (creation) operator; \i denotes the chemical 
potential and h is the external magnetic field. We shall consider h = in much of, and non- 
half-filled bands throughout this paper. Eq. ( |2.1| ) represents the Hamiltonian on N sites 
with open boundary conditions (OBC). For periodic boundary conditions (PBC), extend 
the first sum to N and identify Cn+i = CV 

This model is solved byBethe AnsatzEa for periodic and by reflection Bethe Ansatz for 
open boundary conditions^, and the energies of the ground and excited states are given 
as the solution of certain integral equations. Here we do not repeat these equations from 
the literature but emphasize the main physical results on passing from periodic to open 
boundary conditions. The basic idea of the reflection Bethe Ansatz is to superpose, in order 
to obtain the solution for open boundaries, the Bethe solution of the periodic system with its 
reflection at the originS. In a similar way, one can pass from a particle on a ring to a particle 
in a box in elementary quantum mechanics, from a Luttinger liquid with periodic boundary 
conditions to one with open boundaries^!, and from conformal to boundary conformal field 
theory!. A many-particle eigenstate with N c — + N± electrons out of which N s = iVj have 
spin projection -1/2, is parametrized by the set of rapidities (k, A for charge and down-spin, 
basically a generalization of the wavenumbers of free electron states to include the correct 
scattering phase shifts) of the occupied statesS. While for PBC, rapidites +k and —k give 
linearly independent solutions, for open boundaries only one sign (say +) for these rapidites 
is allowed (as in elementary quantum mechanics), but the spacing of solutions is half as 
big and their density twice as big as for PBC. The equations determining the energies 
are formally different for open and periodic boundary conditions. However, they can be 
transformed into each other (cf. below for the entries of the dressed charge matrix), and 
consequently, their solutions are equal. The energies thus are independent of the boundary 
conditions, as are the velocities of the low-lying collective charge and spin modes (holons 
and spinons). 

The critical properties are determined by the low-energy excitations, more specifically 
by the quantum numbers associated with the operators whose correlations at criticality we 
wish to determine, and a "dressed charge matrix" This 2 x 2-matrix contains the 

effective renormalized coupling constants within and between the low-energy (charge and 
spin) sectors of the Hilbert space of the Hubbard model, and therefore directly determines 
the non-universal critical exponents. It is defined as 

2?(p,b) — ( Zee Z cs \ ^ 
\ Z sc Z ss I 
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The superscripts p, b label periodic and bounded systems. In the absence of a magnetic field, 
jj(p>&) j s completely determined by its first element as 

Z M = ( W J_ j ( 2 - 3 ) 

where k is a cutoff determined by the particle density, again identical for periodic and open 
boundary conditions. 

£( p,fe )(fc ), for a periodic system, obeys the integral equationS 

f<p)(Jfe) = i + J_ [ ° dk'cos(k')K^(smk- smk')^ p \k') (2.4) 

27T 

with the kernel 



g M M= /"^ g^§ cosM. (2.5) 
Jo cosh(cj{//4) 

For open boundary conditions, on the other hand, the integral equation for £^(fco) reads 

(jfe) = i + J_ /" °dA/ cos(A; / )K (fc) (sin k, sin fc')£ {6) (#) (2.6) 
2ti Jo 

with the kernel 

K® (z, z') = (z - z') + (z + z') . (2.7) 

Using K^ b \—z, z') (z, z'), a consequence of ifW(-z) = ifW(^), Eq. ( |2.5| ), we find 

that can be continued to negative with ^ b \—k) = ^ b \k). We can then change 

variables k — > — in the contribution of the second term of to the integral in ( |2.6| ), and 
find 

^(6)^) = £(p)^ and thus Z (6) = z (p) _ (2.8) 

Despite the differences between Eqs. ( |2.4|) , (|2.5|) , and (|2.6|), ( |2.7| ), the dressed charge matrices 
are identical for open and periodic boundary conditions. This translates the fact that they 
are a property of the Hamiltonian only and independent of boundary conditions. The 
boundary effects are contained completely in the representation of the boundary critical 
exponents in terms of the entries of the dressed charge matrices. This representation depends 
on the boundary conditions. 

With a magnetic field, one has a set of four coupled equations for the four entries of 
Z^P' b \ containing different kernels. However, all symmetries used, and all transformations 
carried out above continue to be applicable. Consequently, the dressed charge matrices of 
open and periodic systems are identical also in a finite magnetic field. 

The rapidities k, A, describing a general quantum state of the ID Hubbard model depend 
on two sets of quantum numbers I c j and J s j of integers or half-odd integers, parametrizing 
the solutions (the "occupied states"). There are certain parity rules for the I c j, I St j, de- 
pending on N c , N s being even or odd. In the ground state, their distributions are filled up 
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to a cutoff, a pseudo-Fermi number. Acting with an operator on a state of the system will 
change the distribution of rapidities. Both for the charge and spin channels, in the periodic 
systems, there are three types of excitations: particle addition (often also termed "charge" 
excitations!!), current excitations, and particle-hole excitations. The single-particle opera- 
tor Cl s , e.g., adds a particle, i.e. AN C = 1, and, depending on spin, AN S = for s =| 
or AN S = 1 for s =1. Changing N C}S by unity changes the I c j and/or I s j between inte- 
gers and half-odd integers. This backflow is accounted for by current quantum numbers 
D c = (AN C + AN s )/2 mod 1, and D s = AN s /2 mod 1. The role of these excitations, and 
of the particle-hole excitations, can best be seen by considering, e.g. the density operator 
J2s c n, s c n,s = N~ l J2k,q,s c k+q,s c k,s- This operator does not change the particle number, i.e. 
AN C = AN S = 0. For small q, this operator creates particle-hole excitations, i.e. changes 
the distribution of the positive (or negative) I c j or I s j by some AI C or A/ s . If q is not small 
and rather a multiple of 2k p, particles are transferred across the Fermi surface which will 
generate a persistent current in the system. These large-g excitations therefore change the 
symmetry of the distribution of the J CJ - and I s j, and therefore have a finite D c s . Permissible 
values of the D C)S follow from the expressions above. The finite D C)S values then generate 
the 2/cf, etc. excitations. 

This set of quantum numbers is characteristic for an operator 0(x). Conformal field 
theory then determines its correlation functions as!§ 



a(D c , D s ) exp(—2iD c k F ^x) exp(— 2i[D c + D s ]k F ^ 



The anomalous dimensions of the operator 

. i If Z SS AN C - Z CC AN S \ 2 A , . 

A c=g(± 2detz - + Z cc D c + Z sc D s j +AJ C ± (2.10) 

± If Z SC AN C - Z CC AN S n \ 2 . 

Af = - =F TTTTT? + Z CS D C + Z SS D S ) +AIf (2.11) 



2 V 2 det Z 

are determined by the entries of the dressed charge matrix and the changes in the quantum 
numbers AiV C)S and _D CjS which the action of the operator generates. 

For bounded systems, only positive rapidities are relevant, cf. Eq. ( |2.6| ), there is only one 
pseudo-Fermi point for each channel, and there are no current excitations. Consequently, 
there are only charge excitations AiV c , AiV s , and particle-hole excitations AJ C , A/ s . The 
backflow terms, resp. asymmetries of the rapidity distributions, D c , D s are absent. 

The complete, explicit behavior of correlation functions is rather complicated close to 
an open boundary, cf. Eq. (8) of Wang et a/.0, for an example for a single channel (charge 
or spin). Asymptotically, however, the correlation functions Go of an operator O simplify 
considerably and behave, close to an open boundary at x = 0, 

1. in time (for x±,X2 ~ 0, and vt ^> Xi,^) as 

G (x l ,x 2 ,t) ~r 2x °> h (2.12) 
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2. in space (specifically: x±,t ~ 0, x 2 Xi,vt) as 

G (x 1 ,x 2 ,t)^x^ Xo ' b+do) . (2.13) 

(Only the conformally invariant, i.e. slowly varying part of the correlation function is 
given. This must be multiplied by a factor oscillating with an appropriate multiple of 
k F x 2 .) 

3. in temperature (xi, x 2 ~ 0, uj <C T) 

Go(xi,x 2 ,u,T)^T 2x o^ 1 , (2.14) 

where the Fourier-transformed (in time) correlation function is used. That the bound- 
ary critical exponent xo,b, describing temporal correlations, comes in for temperature is 
most easily rationalized from the perspective of the Matsubara formalism of imaginary 
times (frequencies) in many-body physics. 

The exponents xo,b arid do are the boundary critical exponent and bulk scaling dimension, 
respectively. The different exponents for spatial and time decays, and the appearence of 
a new boundary critical exponent describing the time correlations, translate the combined 
influences of the interactions and of broken translational invariance, on the low-energy prop- 
erties of interacting ID electron systems close to boundaries. In terms of the dressed charge 
matrix Z, and the Quantum numbers AN C , AN S , AI C , and AI S associated with the operator 
O, xo,b is given hy^ 

x ,b = ^- X AN) T • (Z- ! AN) + AI C + AI S (2.15) 

with (AN) T = (AN C , AN S ), and the superscript T denotes the transpose. The bulk critical 
exponent is 

d = A+ + A c - + A+ + A;. (2.16) 

Finally notice that when the system has boundaries on both sides, i.e. is of finite length L, 
the power laws of space and time are changed to power laws of sin(7rx/2L), resp. sm(irvt/2L) 
with the same exponents. 



B. Luttinger liquid theory 

At low energies, models of correlated ID electrons such as the ID Hubbard model, 
reduce to Luttinger liquidsHi. Haldane's Luttinger liquid conjecture! postulates that, in a 
low-energy subspace, gapless one- dimensional quantum systems are described by an effective 
Luttinger HamiltonianEj 

#=^E Jdx{v»Krc 2 nl(x) + ^[^^) \ . (2.17) 
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$!, and 11^ are phase fields and their conjugate momenta for charge [y = p) and spin [y = a) 
fluctuations 



and 



1 dQJx) 



in e ~a\p\/2-ipx 

— r X] K(p) + u -(p)\ 



(2.18) 



7T (9x 



777- , . 

n ^ x) = ^rr^rz with Qu(x) = t -Y J 



-a\p\/2—ipx 
P 



[v+(p) -v_{p)\ 



(2.19) 



which obey bosonic commutator relations, as do the long-wavelength fluctuation operators 
v r (p)- The K v are the renormalized coupling constants in the charge and spin sectors. 
Zero-modes (p = 0) do not influence the dynamics, and have been neglected. 

Correlation functions and their critical exponents can then be calculated by 
bosonizationli 



Arkpx 



^ rs (x) ~ lim — == exp — = [r$ p {x) - Q p (x) + s {r<f> a (x) - & a {x)}\ 



^2.20) 



The chirality label r = ± here and in Eqs. ( [2.18 ) and ( 2.19| ), refers to right- and left-moving 
fermions, with k m rkp- Consider a general local operator 



{m} (x) = O m+vm+l>m _^ m _ L (x) = ^(x)^(x)*™-\x)^(x) 



(2.21) 



where positive (negative) exponents label powers of creation (annihilation) operators, and 
higher powers of the operators are understood to be point-split, e.g., (m+f > 1) 



TO +T 



(2.22) 



i=l 



The two-point correlation function of such an operator decays as 

exp(— ikpxJp) 



(O {m} (xt)O\ m} (00)) 



(x — Vpt) 2d P (x + V p t) 2d P (x — V a t) 2d ° {x + V a t) 



\2d- 



[2.23) 



where the scaling dimensions 

1 r 



2d+ 



2d~ p 
2dZ 



1 r 



(AN P + J P Y + (K p - 1)J P + (K; L - 1)AN P 
(AN P - J p f + (K p - 1) J p + (K; 1 - 1)AN P 



[2.24) 



1 r 



1 r 



(AN a + J a f + (K ff - l)J a + (K- 1 - l)AN a 
(AN a - J a f + {K a - 1 ) J (j + {K~ l - l)AN a 
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are determined by the two Luttinger liquid parameters K u and the number of charge and 
current excitations in the charge and spin sectors, created by the operator 0{ m }(xt) 

AN p = ^ m r,s , Jp = Yl rm r,s , AA^ = Y Sm r,s , J<r = $Z rSlll r,s ■ ( 2 -25) 
r,s r,s r,s r,s 

In addition, there are particle-hole excitations with chirality r in each sector, given by 

Pr {x) = J2 #Ux)* r .(x) , a r (x) = £ s^ rs (x)% s (x) . (2.26) 

s s 

If an operator O differs from O above by AI r ^ powers of u r (x)-type particle- hole excitations, 
the exponent of its correlation function will be increased as 

d r v -^cT u + Al r , v . (2.27) 

Up to now, these critical exponents have been derived case by case, by an explicit bosoniza- 
tion calculation. Our Eqs. (|2.23| ) and ( p.24[) show that, in complete analogy with conformal 



field theory, it is possible to give a general construction rule for the correlation functions, only 
based on the knowledge of the K v and the quantum numbers associated with the operators. 

This picture is well-established for periodic systemsi. For open systems, the same basic 
ideas as outlined above, for integrable models and boundary conformal field theory, continue 
to hold: as a consequence of the boundary conditons, right- and left-moving fermion fields 
are no longer independent, 

# +>s (x) = , (2.28) 

and the physical fermions are described as superpositions of a chiral fermion with its reflec- 
tion at the originEHll. The basic bosonization formula ( |2.20| ) is then replaced by 



~ lim ^= exp (-L [$ (x) + sQ+rix)]) , (2.29) 

where the fields $ +>u (x) are obtained from Eqs. Q2.18Q or (|2.19|) by dropping the V-(p)- 
operators. Due to the open boundary conditions, the full correlation functions of these 
operators are quite involved. While the single-particle Green's function has been calculated 
by othersEHi3, general correlation functions have not been published to date. Here, we give 
simple rules how the exponents describing the asymptotic decay of the correlation functions 
of more complicated operators, close to an open boundary, can be constructed in complete 
analogy to the recipes of boundary conformal field theory. 

In the presence of open boundaries, the general operator 0{ m }(x) is represented in terms 
of only one type of chiral fermion, say r = +. There is thus only a single Fermi point, and 
no current excitations can be defined. There are only charge and particle-hole excitations. 
We therefore have — > OfK where the superscripts refer to boundary conditions 

Ofi } (x) = Kr +m -'(*)K! i+m - l (x) = <f T (*)<fV) (2.30) 

with the same conventions as for oj^i, Eq. ( |2.21| ). The first equality directly translates the 
bulk operator, in terms of its bulk quantum numbers, into a boundary operator, while the 
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second equality only uses the quantum numbers defined for the boundary operator: AN S is 
the number of spin-s particles the operator O adds to the system. One now can calculate 
the correlations functions of Of\{x) by the methods developed by othersM. 

Close to an open boundary in a semi-infinite system, the complicated expressions simplify 
considerably in the limits t>i and x ^> t where x is close to the boundary. From these 
limits, a boundary critical exponent xo,b can be defined for each operator 0{ m y(xt), and the 
behavior of the correlation functions as functions of the variables x, t, and T is formally 
identical to Eqs. ( |2.12| )-( p.l4| ). For a Luttinger liquid, the boundary critical exponent xo,b 
then is given by 

o (p) , (p) M _ (AA^) 2 / 9 01 \ 

■Z%0,b — %o,b ' X 0,b > X 0,b ~ ■ [Z.Ol) 

Formally, AN U is calculated in the same way as above, Eqs. (|2.25|) , when the operator is 
represented in terms of the bulk chiral fermions, resp. with the correct number of m +tS 
and m_ )S = when it is set up directly as a boundary operator. Physically, an operator 
which generated current excitations in a periodic system, will now generate particle-hole 
excitations. The spatial decay of the correlations of 0(x) close to a boundary involves the 
bulk scaling dimension do = J2r,^d r u of the operator 0(x) as in fl2.13p , and the d r v are taken 
from Eqs. ( |2.24| ). The temperature variation again is determined by the boundary critical 



exponent xo,b as in Q2.14|) . Bosonization of a Luttinger liquid with open boundaries there- 
fore precisely reproduces the structure of correlation functions, and the boundary critical 
exponents that boundary conformal field theory extracts for integrable systems. 

Notice further that for Luttinger liquids with open boundaries, the value of K p again 
comes out identical to that for periodic systems, translating the fact that it is a property of 
the Hamiltonian, resp. thermodynamics of the system, and as such independent of boundary 
conditions. Besides determining the boundary critical exponents of ID fermions with open 
boundaries, our work therefore provides rules, in terms of K p , how to connect bulk and 
boundary critical exponents of different correlation functions. Below, we shall discuss an 
application of this procedure. 

Both methods, boundary conformal field theory and the Luttinger liquid, must lead to 
identical critical exponents. The relation between the entries of the dressed charge matrix, 
and the Luttinger liquid K p in the absence of magnetic fields (then K a = 1 for spin-rotation 
invariance) is 

eth) = 2K P , (2.32) 

both for bounded and for periodic systems. 

In a finite magnetic field, the situation is far more complicated. In the ID Hubbard 
model, there are (at least) three important effects of a finite magnetic field 

1. All entries of the dressed charge matrix acquire /i-dependent corrections. This implies 
that the Luttinger coupling constant K p , to the extent that it still makes sense, would 
depend on h and therefore that charge-spin separation is violated, even in a low-energy 
subspace. Explicit expressions have been given by Frahm and Korepiniil. Specifically, 
Z cs 7^ now. 
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2. The correction to the matrix element Z ss is particular in that it is logarithmic in h 

z "^7l( 1 + lto(b))' (2 ' 33) 

while all other corrections have weaker /i-dependences. h c is the critical field for a 
completely spin-polarized state. For the ID Hubbard model, it has been calculated by 
Frahm and Korepin0 and varies from a finite constant at[/ = 0toal/[7 behavior as 
U — > oo. For non-interacting ID electrons, one has h c = e(2kp) where kp is the Fermi 
wavevector for h = and e(k) is the single-particle dispersion taken from the bottom 
of the band, e(0) = 0. Thus h c = 2t[l — cos(2kpa)] for ID tight-binding electrons. 

3. There are different Fermi wave vectors for up- and down-spin electrons, kp^ ^ kp^. 

Quite generally, the mapping to a Luttinger liquid is problematic and cannot be carried out 
any more. If one is interested in boundary critical properties of ID integrable systems, the 
conformal field theory approach (|2.12| )- (|2.15|) must be preferred. 

In the (apparently realistic) limit h -C h c , the situation simplifies, however, and a map- 
ping to a Luttinger liquid can still be done to logarithmic accuracy. To this accuracy, the 
only effect of the magnetic field is the correction to Z ss , Eq. (|2.33| ). This translates in a 
/i-dependent deviation from unity, of the Luttinger coupling constant 

K - = { 1 + Ii^)' Ml+ 2j^Jh)- (2 ' 34) 

We obtained this equation by evaluating the finite-size corrections to the ground state 
energyH, 



£(AN,D)-£ = ^ ' 



v p (A+ + A-)+v (A+ + Aj)\ (2.35) 
from states with finite AN S but vanishing AA C = D c = D s = 0, and using ( |2.10|) and (p. 11 



E(AN s ,0,0,0)-E = ^v a (^-^ . (2.36) 

Its second derivative with respect to AN S can then be identified to the inverse susceptibility 
of a Luttinger liquid and gives ( |2.34| ). To logarithmic accuracy in h, K p is unchanged, and 
one can use also the expressions for the (boundary) critical exponents of Luttinger liquids. 
A relation similar to Q2.34p is also known for the ID Heisenberg model0, and a correction to 
the magnetic susceptibility of the ID Hubbard model, logarithmic in h and consistent with 
( |2.34j ), has been derived by Kawano and Takahashiil. 



III. A NEW UNIVERSALITY CLASS: BOUNDED LUTTINGER LIQUIDS 

Universality classes of critical behavior are usually defined in terms of the scaling relations 
between the critical exponents. From the general expressions ( |2.12j )-( |2.14^ ) above, it is clear 



12 



that these scaling relations on the boundary of a Luttinger liquid are different from those 
in the bulk. Luttinger liquids with open boundaries therefore realize a new universality 
class of (boundary) critical behavior which we propose to call bounded Luttinger liquids. 
More specifically, following Haldane's statement for the periodic systems!, we conjecture 
that one- dimensional quantum systems with gapless degrees of freedom and open boundaries 
form bounded Luttinger liquids. They comprise interacting electron systems, of which we 
discuss an example here, but also ID bosons or spin chains£2l. Bounded Luttinger liquids are 
boundary conformal field theories with central charge c = 10. Their nonuniversal exponents 
depend on one open coupling constant per degree of freedom K v which is the same as in the 
bulk. 

Here, we evaluate explicitly a variety of important boundary correlation functions. This 
will display the new set of scaling relations characterizing the bounded Luttinger liquids, 
and thus back our claim of a new universality class. 

The single-electron Green's function (spin-f ) is 

v ' ' 1 M mM ;n ;/ \ sm(k F x)x x *> b d * \ sm(k F x)x~ l ~ a *- b v ; 

where the last equality transforms our notation to the familiar Luttinger liquid notation 
with Green's function exponents a. It involves the following excitations (both in terms of 
the Bethe ansatz for the Hubbard model, and in terms of a Luttinger liquid) 

AiV c = 1, AN S = 0, AI C = AI s = (3.2) 
AiV p = 1, AiV ff = 1, AI U = (3.3) 

In the absence of a magnetic field (Z cs = 0), its boundary critical exponent is 

2x » = %0 {3A) 



and 



x*, b + rig = \ { 3 Z l s + Z 2 CC + Z 2 SC + Zl - 2Z CC Z SC ) . (3.5) 



4 V (det Z) 2 
Using (|2.32| ), the Luttinger liquid notation obtains: 



and 



The same expressions are obtained, of course, when using directly Eq. ( [2.31| ). Photoemission 
experiments measure the Fourier transform of Gq,(t), the local spectral function 



13 



p(x, lo) = --lmG(x, uj + fi) ~ \u\ a *> b , (3.8) 

7T 

where fi is the chemical potential. For repulsive interactions K„ < 1, we have 



<4 >b > a %p = -\ K p + — -2j (3.9) 

which suggests that the pseudogap in the local density of states is deepened by the presence 
of boundaries. We will use these expressions in Section |TV] below to discuss photoemission 
experiments in organic conductors. 

Due to its relevance to photoemission experiments, the conformal field theory predicition 
( p.8|) has been checked for the ID Hubbard model. No exact calculation of the entire spectral 
function is available. Conformal field theory therefore provides the only exact prediction of 
properties of the spectral function, though only of its exponents. The question therefore is (i) 
to what extent other methods produce exponents consistent with conformal field theory, and 
(ii) if the deepening of the pseudogap in the local density of states is found, too. These issues 
were addressed recently, using perturbative Hartree and numerical DMRG calculations^. A 
suprising result of the Hartree calculations was an enhancement of the density of states, 
instead of a suppression, close to the chemical potential and close to the open boundary. 
DMRG was not able to calculate the frequency dependence of p(x, u). For the zero- frequency 
matrix elements, it asymptotically verified, however, a decay with system size as L~ a<s, ' b , as 
expected from conformal field theory. Moreover, at finite L, the matrix element approached 
the predicted power-law from below. Both facts suggest that the depth of the pseudogap is 
increased by the presence of open boundary conditions. However, that work also suggests 
a word of caution: for small U, the asymptotic power-laws were approached only on very 
large lattices (L > 100 sites)! implying very small energy scales for (|3.8j ) to be observed. 
How these scales depend on the interaction strength or range, and if the pseudogap is over- 
or underestimated at smaller L (the DMRG data seem to suggest the second possibilityH3) , 
deserves further study. It could also be interesting to perform an exact calculation of the 
entire spectral function at U = oo as has been done successfully for periodic systemfi. 

The boundary critical exponent of the Green's function also determines the temperature 
dependence of the tunneling conductance G through a weak link between two Luttinger 
liquids (LL) or between a Luttinger liquid and a normal metal (FL) 

' ~ i +1 l 1 } (3.10) 

In magnetic fields, and for OBC, the exponents controlling these temperature dependences 
depend in first order on \n(h c /h), unlike periodic systems where the dominant dependence 
of the Green's function exponent is of second order in \n(h c /h). This also holds for the 
related situations where particles scatter off a weak impurity in a periodic ring^l because 
here, the temperature dependence of the conductance is controlled by the 2fcp-part of the 
density correlations, 

G{T) = G [l - const. T K " +K °- 2 ] . (3.11) 
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Here, a finite h will increase the exponent while in the tunneling case, the exponent will 
be decreased. If K p < 1, in both cases the effect of a finite magnetic field will be to 
offset the interaction effects on the exponents contained in K p , and to produce conductance 
variations more similar to free electrons. Precisely this effect has apparently been observed 
in semiconductor quantum wires in magnetic fields^. Notice that in such structures, E F is 
small, typically some meV. Consequently, h c is also small and the magnetic field effect on 
K a could well be big. 

The density-density correlation function 

G n (xt, x 2 , t) = (n(x 1 ,t)n(x 2 , 0)) , (3.12) 

involves the following excitations 

AN C = AN S = 0, AI C = 1, A/ s = or AI C = 0, A/ s = 1 , (3.13) 
AiV p = 0, AN a = 0, AI P = 1, AI a = 0. (3.14) 

Notice that for open boundary conditions, there are no current excitations while in the bulk, 
the density operator can involve current excitations and therefore oscillate with wavevectors 
that are even multiples of kp. Consequently, there is a single time behavior while the spatial 
dependence contains the usual q ~ 0, 2k F , 4k F , . . . components. We have 



Xn,b = 1 (3.15) 



and 



G n (0,0,t) ~r 2 , (3.16) 
G n (x, 0, 0) ~ x- 2 + sin(2k F x)x- 3 - ± ^ R + sm(Ak F x)x- 1 ' 2K " . (3.17) 

In the bounded system, there is now a universal boundary critical exponent x Ut b = 1, while 
the time dependence was non-universal in the periodic system. The nonuniversal decay in 
x-direction, and the differences in the exponents of the various momentum components, are 
a consequence of the different bulk conformal dimensions of the 2/cF-harmonics of the density 
operator. 

Friedel oscillations appear in the system because of the breaking of translational in- 
variance by the boundaries. They have been studied in detail by Bediirftig et aZ.,0. We 
determine their spatial decay aa!2l 

.... . sin(2k F x) . sin(4:k F x) 

(n(x)) = n - A, } - A 2 y / } (3.18) 

X 2 X 

where we have used the Luttinger liquid language^. Notice that they are governed not by 
boundary critical exponents but by the bulk exponents of the 2k F - and 4/c^-parts of the 
density operator. The boundary exponents of the density correlations, are different from 
those of the Friedel oscillations, however. DMRG calculations of the Friedel oscillations 
are in impressive agreement with conformal field theory predictions both for spiniest and 
spin-1/2 fermion systems. 
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The longitudinal and transverse spin-spin correlation functions involve the following ex- 
citations, respectively, 



G z a (xi, x 2 , t) = (S z ( Xl ,t)S z (x 2 , 0)> , (3.19) 

S z (x,t) = -K(x,t) - n;(x,t)], 

AN C = AN S = 0, AI C = 1, AI S = or AI C = 0, AI S = 1, 
AiV p = 0, AN a = 0, A/ p = 0, A/ ff = l;, 
G^(a; 1 ,x 2 ,t) = ^-(x^t)^^^)) , (3.20) 
S + (x,t) = v\(x,t)V l (x,t), 

AN C = 0, AN S = 1, AJ C = AI S = 0, 
AiV p = 0, AA^ ff = 2, A/ p = 0, AI a = . 

In the absence of a magnetic field, they decay with the same exponents as the long- 
wavelength and the 2/cp-components of the density correlations. The enhancement of mag- 
netic over density correlations which we expect in the U > 0-Hubbard model, in the bulk 
is due both to logarithmic corrections and to prefactorslil. A similar enhancement in the 
correlation functions, is expected close to the boundary. Notice, however, that the static 
Friedel oscillations in the density will dominate there, and the importance of the magnetic 
correlations in the system will be limited either to dynamical measurements, or to the bulk 
region. It is also important to notice that, in the Luttinger liquid picture, it is the spe- 
cial value K a = 1, embodying spin- rotation invariance, which makes the boundary critical 
exponents of the longitudinal and transverse spin correlations identical. Formally, we have 

4, 6 = 1 and x^ b = l/K a . (3.21) 

Finally pairing correlations are important. The singlet pairing correlations are 

Gf s (x u x 2 ,t) = (* T (a; 1 ,t)*;(a; 1> t)*[(x 2> 0)*{(ar 2 ,0)) , (3.22) 
AiV c = 2, AN S = 1, AI C = AI S = 0. 

AiV p = 2, AN a = 0, AI V = . (3.23) 
Their boundary critical exponent is 

=m = k (3 - 24) 

so that they decay as 

G ( sl(0,0,t) ~r 2/K ? (3.25) 

G { sl(x, 0, 0) ~ x -^- P -^ (3.26) 

They will dominate for attractive interactions [K p > 1) only. Then, however, the presence 
of a boundary will enhance them over their bulk values. 

Triplet pairing correlations have three components (S z = 1, 0, —1)) involving the follow- 
ing operators and excitations 
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G? s (x u x 2 ,t) = (^ T (x 1 + l,t)^ T (x 1 ,t)^{(a; 2 ,0)^ t T (x 2 + l,0)) , (3.27) 
AiV c = 2, AN S = 0, AI C = AI S = 0, 
AN p = 2, AA^ CT = 2, AI„ = , 

GP s ( Xl ,x 2 ,t) = (^(x^^t^lix^O^x^O)) , (3.28) 
AN C = 2, AiV s = 1, AJ C = A/ s = 0, 
AiV p = 2, AiV CT = 0, AI V = , 
Gfe 1} (x ls x 2 , t) = (^(x! + 1, t)^(xi, t)*I(x 2 , 0)^j(x 2 + 1, 0)) , (3.29) 
AN C = 2, AN S = 2, AI C = AI S = 0. 
AN p = 2, AN a = -2, AI„ = , 



Using the quantum numbers of the excitations, we find for S z = ±1, 

= |4) + 1 = ^ + 1 (3 ' 30) 

This is different from the value of S z = which is equal to that for singlet pairing, and given 
in ( 13.241) . Superficially, this would imply a breaking of spin-rotational invariance which is 



inconsistent with the model Hamiltonian. The solution of the puzzle is contained in Eq. 
( |2.28| ). In a periodic Luttinger liquid, the triplet pairing operators with spin projection 



S z = +1 and are represented as 
O ts ,+i(x) = ^+, T (x)^_, T (x) and O TS , (x) = [*+, T (a;)*_ a (a;) + * +)i (a;)tf_, T (aO] 



(3.31) 



The linear dependence of right- and left-moving fermions, Eq. ( |2.28| ), in the presence of open 



boundaries clearly affects the operator T5 +1 (x), and is accounted for correctly by our rules 
for the critical exponents. On the other hand, the individual contributions ^/ + ^(x)^/-^(x) to 
Ots,o(x), underlying the critical exponent (|3.24j ) would not seem to be affected by ( |2.28|) be- 



cause different spin projections are involved. However, when adding both terms in Ots,o(x), 
it is clear that they add up to zero when using fl2.28j ). Therefore, the pref actor of the term in 



the correlation function G^' s (x\,X2,t) carrying the boundary exponent (|3.24|) must vanish, 
and only the next higher term survives. We therefore conclude that the boundary critical 
exponent of the triplet pairing correlations is 

x TS,b = x Tsl = ToTTT + 1 = TT" + 1 • (3.32) 



p 



The singlet correlations are not affected by this argument and ( |3.24j ) continues to hold. 
Note that xxs,b = xss,b + 1 while their bulk scaling dimensions (and therefore their critical 
behavior in periodic systems) are equal. 

Gapless ID quantum system with open boundaries form a new universality class which 
we propose to call "bounded Luttinger liquids" . As we have shown earlier0, they are c = 1 
conformal field theories, or products thereof, and they posses critical exponents different 
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from bulk Luttinger liquids. With bulk Luttinger liquids, they share the fact that the scaling 
relations between the nonuniversal critical exponents are parametrized by one renormalized 
coupling constant K v per degree of freedom v. These coupling constants are characteristic 
of the underlying Hamiltonian and the same irrespective of the boundary conditions. Unlike 
bulk Luttinger liquids, bounded Luttinger liquids have no current excitations. This gives 
universal boundary critical exponents to correlation functions of operators involving current 
excitations in the bulk, where they had nonuniversal exponents. Examples are the 2kp- 
and 4fc^-components of the density correlations. This should also apply to the 3/cp-Green's 
function which would be another interesting example that might be worked out. Also, as 
we have shown, the linear dependence of right- and left-moving fermions introduced by the 
boundary, may produce different boundary critical exponents for operators which have the 
same bulk scaling dimensions. 

Finally, while charge-spin separation may be present in a bounded Luttinger liquid, it 
will be exceedingly difficult to observe it: all excitations are localized by the boundaries, 
and there are no propagating solutions. 

When expressed in terms of Kn, K a , our single-particle Green's function takes the same 
form as those found elsewheretBcJ. The suggestion that a new universality class is realized 
by correlated fermions with open boundaries, was also made by Mattsson et a/.0. 



IV. APPLICATION: THE PHOTOEMISSION PUZZLE IN THE BECHGAARD 

SALTS 

Here, we discuss one way in which these results can be put to use. The family of organic 
conductors (TMTSF) 2 X (where TMTSF stands for the molecule tetramethyltetraselena- 
fulvalene, and X typically is PF%, CIO4, ReO^, etc., the "Bechgaard salts") is one of the 
candidates for the realization of a Luttinger liquid in quasi-lD electron systems. There is 
general experimental evidence for strong anisotropy of the electronic properties and strong 
electron-electron interactions. More specifically, several experiments have produced evidence 
in favor of Luttinger liquid behavior: 

1. The NMR spin-lattice relaxation is consistent with a Luttinger form Tj~ ~ T + T Kp . 
Values K p ~ 0.15 were sugg estedfl. 

2. Optical conductivity does not directly probe the Luttinger liquid but its leading ir- 
relevant operators which can relax momentum to the lattice. Assuming that these 
stem from electron-electron Umklapp scattering in a quarter-filled band, theory pre- 
dicts <y{o*>) ~ uj 16Kp ~ 5 on the high-frequency tail of a mid-infraredDeakll. Experiments 
observe power-law behavior in this region, leading to K p ~ 0.23E3. 

3. Under a pressure of 9 kbar, the resistivity perpendicular to the chains shows power- 
law behavior in temperature^. Theoretically, the exponent is related to that of the 
single-particle Green's function, giving values K p ~ 0.25 . . . 0.3. 

Although the values of K p suggested by the experiments do not agree precisely, the proper- 
ties probed apparently are consistent with a Luttinger liquid. There is, however, one type of 
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measurement whose results systematically disagree with Luttinger liquid predictions: pho- 
toemission. 

The spectral function of a Luttinger liquid generically exhibits two dispersing signals, 
due to charge-spin separation, whose singularities depend on i^ p ,i'0Jli. For K p ~ 1/3 ... 1/4, 
one expects a divergence dispersing with the charge velocity and an exponent — 1/3. . . — 
7/32 pd — 1/4 while the singularity dispersing with the spin velocity has a smaller exponent 
— 1/6 ... + 1/16. This is not observed^. The result of the experiment on (T 'MT 'S 'F)2ClO '4 
can be summarized briefly: no charge-spin separation, no dispersion, not even low-energy 
peaks. Instead, there is a broad peak at -1 eV whose tail reaches down to zero intensity at 
the chemical potential!!! The leading edge of the spectra is independent of the wavevector 
k of the photoelectrons, and well fitted by a power law, ~ (— u) a , with exponents a > 1 
over a wide range of energies. Moreover, the appearence of the momentum-resolved spectral 
functions is strikingly similar to that of an earlier experiment with low (or no) angular 
resolutionil. 

We now discuss the possibility that this may be due to impurities on the sample surface. 
The argument involves three main steps which we develop in the following: (i) ARPES 
probes intrinsic properties of the surfaces of the Bechgaard salts; (ii) the surface states 
from which the photoelectrons are ejected, are localized; (iii) if these states are described 
as bounded Luttinger liquids, the A^-values suggested by the bulk measurements together 
with Eq. ( |3.8|) give a good description of the photoemission spectra. 

Very similar lineshapes are observed in the insulating compound (TMTTF) 2 PF^ (sele- 
nium substituted by sulfur) which is believed to be a ID Mott insulator. In photoemission, 
the important difference to (TMTSF) 2 X is that a gap of about 100 meV is observed, 
consistent with the values of the charge gap deduced from other experiments. Very new 
experiments on (TMTSF) 2 Re04 which undergoes a metal-insulator transition due to an- 
ion ordering, observe lineshapes similar to ^TMTSF)2ClO± in the metallic phase and to 
(TMTTF)2PFq in the insulating phase, and furthermore can monitor the gap opening as 
a function of temperature. Photoemission therefore is sensitive to spectral changes at the 
Fermi surface, and observes intrinsic properties of the surfaces of the Bechgaard salts. 

The absence of dispersion in the experimental data then implies that the photoelectrons 
are ejected out of localized states. Then, even an angle-resolved photoemission experiment 
will measure the local density of states, p(x,u) = J dkp(k,u), i.e. the momentum- resolved 
spectral function integrated over a wide range of momenta. The surface sensitivity of pho- 
toemission spectroscopy restricts this conclusion to the sample surface only. In the bulk, 
there may well be the propagating charge and spin excitations which are believed to underly 
the experiments discussed above. (Notice, though, that apart the longitudinal DC resisitiv- 
ity not discussed here, many experiments do not require them explicitly to be propagating.) 
The origin of the localization at the surface is not clear, although impurities, perhaps intro- 
duced by cleaving the samples, and/or the enhanced one-dimensionality, provide plausible 
causes. 

If the bulk is assumed to be a Luttinger liquid, and if we approximate the impurities at the 
surface by open boundaries!!, then the chains at the surface consist of finite segments which 
must be described as bounded Luttinger liquids with the same K p as the bulk. Provided 
the conditions underlying Eqs. (|2.12 ) and ( |3.8|) are fulfilled, the local spectral function is 
obtained as 
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p(x, uj) ~ \uj 
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(4.1) 






*p = 


1/4 





The result for K p = 1/2 has been obtained earlier by Eggert et alS3 } who already speculated 
on the possible influence of surface impurities. In particular for values K p pa 1/3... 1/4 
suggested by optics and DC-resistivity, the low-frequency photoemission weight is well de- 
scribed by fl4.1|) for a bounded Luttinger liquid. Notice moreover that, with photoelectrons 
emanating from localized states, the independence of the results of angular resolution finds 
a natural explanation. 

The conditions on which Eq. ( |4.1| ) is based, imply that most of the ARPES intensity 
must come from states which are dominated by boundary effects. This energy range is widest 
when the impurities are rather closely spaced. In such a case, the discretization of the energy 
levels may no longer be negligible. However, these discreteness effects are likely washed out 
by averaging over many different segments, and one returns to the effective semi-infinite 
bounded Luttinger liquid underlying Eq. (|2.12|) , as demonstrated by Mattsson et alM in 



the related problem of one small segment studied with a poor experimental resolution. A 
more detailed investigation of the spectral properties of a ID correlated electron system with 
random impurities would, however, be useful. 

We therefore come up with a picture of the Bechgaard salts as an inhomogeneous system 
where the bulk would be a standard Luttinger liquid with propagating charge and spin 
excitations, and impurities localizing these excitations at the surface. 



V. OPEN QUESTIONS 

t As a consequence of broken translational invariance, a new set of boundary critical 
exponents appear in correlated ID electron systems close to open boundaries. They define 
a new universality class which we proposed to call bounded Luttinger liquids. In this paper, 
we explored these new scaling relations both by applying boundary conformal field theory 
to the ID Hubbard model, and by using Luttinger liquid theory. Both methods can be 
implemented in a very similar way and, of course, lead to identical conclusions. 

Our study leaves open a number of questions, however. We only investigated open 
boundary conditions. There is the possibility that the boundary critical exponents depend 
on type of boundary, e.g. additional boundary fields (magnetic or chemical potential, etc.), 
boundary operators (e.g. spins, relevant for the Kondo problem@), superconductors, etc., 
and that new scaling relations are associated with them. Such a possibility is suggested by 
a study of a Luttinger liquid coupled to a superconductor^ which finds a local density of 
states p(x,lu) ~ |co>|( A ' p ~ 1 )/ 2 , instead of our ( |3.S| ) and fl3.6|). One could then be left with a 
whole manifold of universality classes depending on the specific type of boundary conditions. 
The dressed charge matrices, and the Luttinger liquid coupling constants K v will, however, 
be independent of these details. 

The bulk - boundary crossover has not been studied systematically yet. While this 
crossover has been studied for the local single-particle density of states in a bounded Lut- 
tinger liquid by Mattsson et a/.0, the important question of how dispersion in a spectral 
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function appears when the length of a bounded electron system is varied from microscopic 
to macroscopic scales, has been left open. In the same way, one would like to understand, 
in terms of dispersion and critical exponents, how the open boundary condition effectively 
emerges in a correlated ID electron system with an impurity. The extension to the spectral 
properties of systems with many random impurities would also be important. 

We used the theory of a bounded Luttinger liquid to discuss the idea that impurities at 
the surface of organic conductors may lead to localized electronic states which would domi- 
nate the photoemission properties of such materialsil. Such an analysis, using the A" p -values 
suggested by experiments probing the propagating bulk excitations, gives a local density of 
states consistent with experiments. Of course, the picture of an inhomogeneous material, 
with localized states at the surface and propagating states in the bulk, is somewhat specu- 
lative. It would predict, however, a definite dependence of the results of an experiment on 
its probing depth. Such dependences should be searched for, and only if they are found, this 
idea should be taken serious. Finally, theoretical information on spectral weight distribution 
beyond the predictions of conformal field theory^, and a clarification of the energy scales 
on which they are accurate, would be very important - even for the periodic systems. 
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